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ABSTRACT 


The pitfalls of the present method of reducing surface pressure to sea 
level are reviewed, and an alternative, adjusted pressure, P, is proposed. P 
is obtained from solution of a Poisson equation over a continental region, 
using the simplest boundary condition along the perimeter or coastline where 
P equals the sea level pressure. The use of P would avoid the empiricisms and 
disadvantages of pressure reduction to sea level, and would produce surface 
pressure charts which depict the true geostrophic wind at the surface. 



1 . Introduction 

The need for standardizing surface pressures measured at elevated land 
stations, so as to essentially eliminate purely hydrostatic variatons, is well 
known (Hewson and Longley, 1944; Saucier, 1955; Wallace and Hobbs, 1977). If 
some such procedure were not performed, the surface pressure map would merely 
resanole en inverted topographic map, due to the approximate 12 mb decrease in 
pressure per 100 m increase in elevation. Synoptic weather influences would be 
largely masked. 

The method presently employed in the US to standardize surface pressures can 
be simmarized as follows (Saucier, 1955; see also Reichelderfer, 1963, Chap. 7 
and Appendix 7.2) If p^ is the sea-level reduced pressure to be plotted, and p^ 
is the actual pressure measured at station elevation above mean level, the 
starting point for p^ is 

Po ' Ps ^ 0.21(mb)-Zg(km)*(T* - T*^)(C) (1) 

where g is gravity, the gas constant for dry air, T* a mean virtual tempera- 
ture, and T*^ an annual normal value of T*. The first term on the right of (1) 
is the hydrostatically correct sea-level pressure that would obtain, for z^ not 
too large, if there were an air column extending below the surface down to sea 
level with a mean virtual temperature T*. The second term is the "plateau" cor- 
rection designed to alleviate inconsistencies noted, on a climatologic basis, 
if only the first term is used. 

In (1), T* is, in its simplest specification, given by 

T* = T*J + (2) 

where T* is the average of the present and 12-hr previous virtual temperature 
at shelter height, and y is one-half the dry adiabatic lapse rate g/Cp, where 
Cp is the specific heat at constant pressure. 
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In seeking a rational scheme for surface-pressure standardization, however, 
this method leaves very much to be desired: 

1) No air colimm exists from z^ down to sea level in general; the quantity 
T* is fictitious. 

2) Above an elevation of 305m, y is taken to be an empirical function of 
surface temperature, involving a diagram for its evaluation. Moreover, this 
variation in y is made a subjective function of geographic location (e.g.. 
Pacific-slope, middle plateau or eastern-slope regions). 

3) The definition of in (2) filters out the diurnal influence of boundary- 
layer warming and cooling upon the reduced pressure, as pointed out by Sangster 
(1960). In many applications (e.g., diagnosis of upslope winds or sea breezes) 

it is necessary to retain the diurnal influence. 

4) In the plateau correction term, z^ is not in all cases the true elevation. 
For a station whose elevation differs "greatly” from that of surrounding stations, 
z^ is the average elevation of the neighboring stations. 

5) To determine precisely the ingredients by which personnel at any par- 
ticular elevated station reduce their own surface pressure to sea level requires 
additional information to be found cnly in non-standard literature references 
which are generally unavailable on short notice. 

6) During conditions of strong temperature contrast over elevated or 
sloping terrain, gradients of sea-level reduced pressure strongly misrepresent 
the actual horizontal pressure gradient at the surface. The reduced isobars 
incorrectly resemble the isotherms then (Sangster, 1960). 

7) Despite the complexity of the method presently in use, discrepancies in 
reduced sea-level pressure for neighboring high stations can still be as large 
as 10 mb in some situations, according to Saucier. 
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An alternative to sea-level reduction of pressure was proposed by Sangster 
(1960); in place of the surface pressure chart one would utilize the geostrophic 
stream function at the surface, and perhaps also the geostrophic potential 
function. The contour shapes and spacings of the former closely resemble 
those of p^ when is small, and otherwise the disadvantages of 1) through 7) 
would be eliminated. However, Sangster's proposal never became operational, 
perhaps because the traditional use of pressure, with its familiar units, would 
have to be abandoned, even over the oceans where there is no problem. 

The present motivations for seeking an alternative to sea-level reduction 
of surface pressure duplicate those of Sangster. However, the alternative pro- 
posed here produces the end result desired by the early architects of the sea- 
level reduction method: standardized pressures whose horizontal gradients yield 
the true horizontal pressure gradients which depict the geostrophic wind at the. 
surface. 


2. The proposed alternative 

The two components of the geostrophic wind, V_, if evaluated at the surface, 

~g 

are 


''g " 


(3) 


where p is density, f is the Coriolis parameter, and subscript s refers to 
evaluation at or very near the surface. Note that, with the horizontal 
gradient operator. 


<V>s * 


'hP. 


over variable terrain, so that a map of p^ alone would not permit to be 

s — 9 

obtained. That is, (v^p)^ is v^p evaluated at the surface, whereas v^p^ is 
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the horizontal gradient of the surface pressure. Therefore we define a horizon- 
tally adjusted pressure, P(x,y,), by 

vP = (Vj^p)g (4) 

Using the hydrostatic assumption, the well known relation (Hess, 1959, Eq. 
(12.9)) between {^^P )5 in (4) and the directly observable quantity vp^ Is 

vP = (Vj^p)g = VPg + Ps^''^s 

Eq. (5) may be re-expressed as 

vP s (v^p)g = 7(Pj + p^gZg) - ZgV(pgg) (6) 

or, using (4) and (6), 

viP - (Pg + Ps92g)] = -ZgV(pgg) (7) 

The divergence of (7) produces the mathematical equation from which P is 
to be solved: 

v2$ s v2[p - (p^ + PgOZg)] = -7*[z^v(p^g)] (8) 

Upon denoting the solution to this Poisson equation by S(x,y), (8) becomes 

= -v*[z^v(pjg)] (9) 

with the result 

P = Pg + Ps92s * S (10) 

Over the oceans, z^ = 0 (since pressure observed at ships' bridges Is re- 
duced to sea level). Then (9) or (7) Indicates that S = 0 upon choosing the cons- 
tants of Integration to be zero. Thus, (9) can be solved more simply by standard 
techniques over son« region whose borders lie at or close to sea level so that 


the Dirichlet boundary condition 

S = 0 (at perimeter and over oceans) (11) 

applies. That is, where i = 0 (10) indicates that for S * 0 we have P = p_ 

5 S 

as desired. This property is missing from Sangster's proposal -- that a 
standardized surface pressure should retain the dimensions of pressure and 
become synonomous with observed pressure over the oceans. 

Having utilized (6) instead of (5), it may be noticed that S will often 
be a relatively small addition to pj9Zg which resembles a sea-level reduction 
term (but using instead of a somewhat larger "sub-surface” density). In the 
next section we show that for a horizontally homogeneous atmosphere S is posi- 
tive (in the neighborhood of 6 mb for z^ = 1km). Only for these reasons does 
the fictitious procedure of reducing pressure to sea level work as well as it 
does. 

An alternative to boundary-condition (11) comes from (7): 

vS = -ZjV(pgg) (12) 

which is a gradient condition. Use of (12) permits the technique to be used 
over any section of terrain, even mesoscale, provided P (or 5 ) is known at one 
point, at least, along the perimeter or in the interior. 

3. Solutions for a horizontally homogeneous atmosphere 
In this particular case p^ = PgCz^) so that 

VPg = ' (3p/3z)vzj . (13) 

Then (7) becomes 

V[P - (Pg + -g(»p/3z)7(Zj2/2) . 


(H) 



Upon approximating 3p/3z with Its average value. Indicated by overbar, over the 
terrain heights Involved, (14) becomes 

v[P - (Pj + Pg gZg) + g(3p/3z)Zg^/2] = 0 

which yields 

P = Pj + PjPZj + g(F/T*)(g/Rj - 7 ) 2 j ^/2 (15) 

upon replacing the density gradient with Its hydrostatic value for a mean lapse 
rate, 7 » and upon noting that the constant of Integration Is zero. Comparison 
of (15) with (10) Indicates that the third tenn on the right of (15), which Is 
positive, represents S. In this special case no further Integration In the 
horizontal Is required to obtain P. 

An expansion of (1) for this specie! case can easily be shown to give es- 
sentially the same result as (15), if the plateau correction term Is omitted. 

Only for this case, then, does the present procedure of reduction to sea level 
make physical sense. In the actual atnosphere, and especially when horizontal 
temperature gradients and terrain heights are large, a property adjusted sur- 
face pressure requires solution of an elliptic (Poisson) equation. Surface 
pressure at any given point cannot. In general, be properly standardized Inde- 
pendently of the surface pressures (and p , z values) at neighboring points. 

4. Discussion 

It might be wondered how It Is possible to adjust surface pressures every- 
where over terrain of whatever Irregularity, so that the derivative of P at any 
point yields the true horizontal pressure gradient, and also maintain the correct 
boundary condition P ® P- “ P^ along the perimeter at sea level (S = 0). The 
answer to this question becomes one of showing that the area average of (9) 

Is compatible with boundary condition (11). From the divergence theorem In two 
dimensions, the area Integral of v*vS on the left of (9) Is the line Integral 



of 7S along the perimeter. The latter Is zero since S Is zero there. The area 
average of the right-hand side of (9) Is zero by the same theorem, since the line 
Integral of z v(p g) Is zero along the perimeter. Hence, compatibility is 
assured, and the average of the forcing function for v^S Is zero. 

The solution for S and P through use of second-order finite differences Is 
not without error, however. To examine Its extent, tests were made with a hypo- 
thetical continent having square sides of length 4400 m and a Nkm plateau In 
the central 15% of the region. The grid Interval was 200 km. A horizontally 
homogeneous atmosphere was prescribed for this and lesser heights, so that the 
true solution was known to be merely the pressure field Imposed hydrostatically 
from greater heights. When a uniform potential temperature, 0, was prescribed 
for the lowest kilometer, the maximum error In P, utilizing (9)-(ll), was 0.06 mb 
while that In p^^ was 2.6 mb (due perhaps to an 1napprop*“1ate choice (1)) 

When a vertical discontinuity In e of 12C was prescribed to occur at z “ 500 m 
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along the maximum slope of the plateau (1.6 x 10 ), the error In P Increased to 
0.46 mb. Since this Is a rather severe test case, the tentative conclusion Is 
that errors of solution for P, apart from uncertainties In Interpolating the 
station network data onto the computational grid, will not usually exceed H mb 
and win be up to an order of magnitude smaller than the errors In p^. An 
additional error of up to a few hundredths of a millibar in P was noted to occur 
at sea level along the coast adjacent to the sloping terrain, when the P * Pq 
boundary condition was applied swnewhat farther out to sea. 

To obtain the surface geostrophic wind requires knowledge of the surface 
density. Even If p^ Is not available to the user, can be estimated to within 
H% accuracy knowing P, T* and z., through use of (10), Ignoring S and utilizing 
the equation of state. One then finds 


P 


/(^*^) 


for use In (3), which becomes 


Ug *• -(pj.")“^3P/3y 
Vg “ (p5f)"^3P/ax 
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An interesting question arises whether on the right-hand side of (9) 

should be smoothed before solving for S. This procedure was not done in the 
example of Sec. 4, and does not seem to be necessary since the solution of a 
Poisson equation is much smoother than Its forcing function. However, Sangster 
did apply such smoothing. 

Another question that arises is how small a height above the surface 
subscript s should refer to in all the preceeding equations. Instrument 
shelter height would seem to be satisfactory, through a somewhat greater height 
within the surface layer might be preferable for better representativeness. 

The proposed alternative to sea-level reduction In surface pressure could 
be implemented much more easily now by the National Weather Service than just a 
few years ago. Thei. recent development of objectively analyzed pressures on 
sane of the 3-hourly surface pressure charts would permit the inclusion of a 
solution of the Poisson equation for P with minimal further effort for an area 
encompassing North America. Presumably, it would be Impractical at present to 
substitute P for p^ In the one-hourly weather reports, since P cannot be obtained 
accurately without the foregoing analysis. 

The use of P Instead of p^ might well encourage comparisons between numeri- 
cal forecasts of surface geostrophic winds and their observations. At present 
there Is little incentive to make such comparisons over elevated terrain, using 
the sea-level reduced pressures for verification purposes, since the latter con- 
tain fictitious elements and since the numerical nx)del may not use the same al- 
gorithm to reduce Its calculated pressures to sea level as do the Individual 
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reportlng stations. 

5. Summary 

A horizontally adjusted surface pressure, P, has been proposed to replace 
sea-levil reduced pressures (p^^) presently In use on the 3-hourly synoptic 
charts. Unlike p^, P has a firm, non-empirical physical basis and is uniquely 
defined. The contours of P show, at a glance, the true direction of the 
geostrophic wind at the surface even over highly elevated terrain. 

Achnowledgment 

This study was supported in part by research grant NAG 5-51 of the 
National Aeronautics and Space Administration. Suggestions by D. Barber 
are appreciated. 


- 10 - 


» N. - 


KEFERENCES 

Hess, S. L, 1959: Introduction to Theoretical Meteorology . Holt, Tinehart and 
Winston, New York, 362 pp. 

Hewson, E. W., and R. W. Longley, 1944: Meteorology Theoretical and Applied . 
John Wiley and Sons, Inc., New York, 468 pp. 

Relchelderfer, F. W., 1963: Manual of Barometry (WBAN) . Vol. 1, U.S. Weather 
Bureau (Dept, of Cotimerce), U.S. Government Printing Office, Washington, D.C. 

Sangster, W. E., 1960: A method of representing the horizontal pressure force 
without reduction of station pressures to sea level. J. Meteor ., 17, 166-176. 

Saucier, W. J., 1955: Principles of Meteorological Analysis . University of 
Chicago Press, 438 pp. 

Wallace, J. M., and P. V. Hobbs, 1977: Atnwspherlc Science . Academic Press, 


New York, 467 pp. 




